First principles study of the solubility of Zr in Al 
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The experimental solubility limit of Zr in Al is well-known. AlsZr has a stable structure DO23 
and a metastable one LI2. Consequently there is a metastable solubility limit for which only few 
experimental data are available. The purpose of this study is to obtain by ab-initio calculations the 
solubility limit of Zr in Al for the stable as well as the metastable phase diagrams. The formation 
energies of several ordered compounds Alj^Zrji.^,), all based on an fee underlying lattice, were 
calculated using the FP-LMTO (Full Potential Linear Muffin Tin Orbital) method. Taking into 
account all the relaxations allowed by the symmetry, we found the DO23 structure to be the stable 
one for AlaZr. This set of results was then used with the cluster expansion in order to fit a generalized 
Ising model through the inverse method of Connolly- Williams. Different ways to consider volume 
relaxations were examined. This allowed us to calculate in the Bragg- Williams approximation the 
configurational free energy at finite temperature. According to the previous FP-LMTO calculations 
the free energy due to electronic excitations can be neglected. For the vibrational free energy of 
ordered structures we compared results obtained from a calculation of the elastic constants used 
with the Debye model and results obtained from a calculation of the phonon spectrum. All these 
different steps lead to a calculation of the solubility limit of Zr in Al which is found to be lower 
than the experimental one. The solubility limit in the metastable phase diagram is calculated in 
the same way and can thus be compared to the stable one. 



I. INTRODUCTION 

The development of methods based on the density 
functional theory^''^ and of the computer power has al- 
lowed to conceive calculations of phase diagram from first 
principles'"^ as an alternative to laboratory experimenta- 
tions. Traditionally, only substitutional effects were con- 
sidered, which was good enough to reproduce the topol- 
ogy of most phase diagrams. So as to obtain a more 
quantitative agreement with experimental data, it has 
been shown more recently that electronic excitationSii as 
well as lattice vibrationaSii could play an important part 
in the relative stability of different phases. We chose to 
test the ability of these first principles methods to predict 
the solubility limit of Zr in an aluminium solid solution 
this part of the Al-Zr phase diagram being interesting 
because of the presence of a metastable phase. 

The Al-richest intermediate phase of the Al-Zr phase 
diagram^ is AlaZr. This compound has the DO23 struc- 
ture which is body-centered tetragonal with eight atoms 
per unit cell, some of these atoms being allowed by sym- 
metry to move along the main axis of the unit cell (Fig. 
CJ. It is stable up to 1580±10°C. 

The solubility limit of Zr in Al (fee) is really low, the 
maximum solubility being 0.083 at.% at the peritectic 
reaction, Liquid-l-ZrAla < — > (Al). The solubilities of 
Zr in both liquid and solid were definitively determined 
by Fink and WillejiS and the assessed phase diagram is 
based on their data. The solid solubility was determined 
from resistivity data and checked by metallography. Solid 
solubilities were also reported by Glazov et alP , Drits 
et ali^ (solubilities determined by means of microstruc- 



tural analysis and electric resistivity measurements), and 
Kuznetsov et ali^ (determined from resistivity, micro- 
hardness and lattice constant measurements as well as 
metallography). The solubilities reported in these last 
studies are higher than the ones of the assessed phase 
diagram. 




FIG. 1: Definition of the structures DO23 (left) and LI2 
(right). 

Supersaturated solid solution of Zr in (Al) containing 
as much as 3 at.% Zr can be prepared by rapid solid- 
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ification. A coherent metastable phase AlsZr precipi- 
tates from the supersaturated solutionis^. This phase has 
the structure LI2 which is simple cubic with 4 atoms 
per unit cell (Fig. 0. This metastable phase can also 
form from the melt as a primary phase during rapid 
solidificationiii^: AlaZr acts as nuclei for solidification 
of (Al) and Zr can thus work as a grain refiner of Al. This 
metastable phase is also responsible for the effectiveness 
of Zr to control recrystallization in Al-based alloys: it 
leads to a more uniform distribution of fine precipitates 
that pins grains and subgrains boundaries. Moreover, 
this phase is really stable against coarsening and redisso- 
lution, all this making it highly desirable. As few experi- 
mental data are available for this phase, it is hard to fit a 
thermodynamic model for it. In such a case, a first prin- 
ciples calculation of the phase diagram should allow us to 
obtain properties that are not available experimentally. 

In order to assess the metastable phase diagram, we 
studied with the same tools and approximations the 
LI2 and DO23 phases: this allows us to check first the 
agreement between the stable phase diagram obtained 
and experimental data, and then to compare it to the 
metastable phase diagram. 

In a first part, we study the stability of ordered com- 
pounds based on an fee underlying lattice in the Al-Zr 
system. The energies of different structures were calcu- 
lated using an ab-initio method, the fuU-potential-linear- 
mufhn-tin-orbital (FP-LMTO). The equilibrium param- 
eters, like the volume, shape of the unit cell, or positions 
of atoms, were obtained. For stable structures they can 
be compared to experimental data. 

Using this whole set of calculations we applied the clus- 
ter expansion to deduce the energy of any structure based 
on the same underlying lattice in the Al-Zr system, ex- 
amining carefully the way to include volume relaxations. 

At finite temperature, the electronic excitations, the 
vibrational free energy, and the configurational entropy 
have to be taken into account. At the end of this part, 
we are able to obtain a thermodynamic model written in 
the same way as in a Calphad approach and to calculate 
the corresponding solubility limits. 



II. GROUND STATES OF ORDERED 
COMPOUNDS 

Formation energies were calculated at absolute zero 
temperature for 26 compounds in the Al-Zr binary sys- 
tem, all based on an fee lattice. Calculations were carried 
out using a full-potential linear-mufhn-tin-orbital (FP- 
LMTO) methodiSiiSiii in the version developed by Meth- 
fessel and Van Schilfgaarde-'-^. The basis used contained 
22 energy independent mufhn-tin-orbitals (MTO) per Al 
and Zr site: three k values for the orbitals s and p and 
two K values for the orbitals d where the corresponding 
kinetic energies were = 0.01 Ry (spd), 1.0 Ry (spd), 
and 2.3 Ry (sp). A second panel with a basis composed of 
22 energy independent MTO with the same kinetic ener- 



gies was used to make a correct treatment of the 4p semi- 
core states of Zr. The same uniform mesh of points was 
used to make the integrations in the Brillouin zone for va- 
lence and semicore states. The number of divisions along 
reciprocal vectors was chosen big enough to ensure a con- 
vergence of the total energy of the order of 0.1 niRy/atom 
for each structure. The radii of the muffin-tin spheres 
were chosen to have a compactness of 47.6% for Al sites 
and 58.4% for Zr sites. Inside the muffin-tin spheres, the 
potential is expanded in spherical harmonics up to Z = 6 
and in the interstitial region spherical Hankel functions 
of kinetic energies — 1 Ry and 3.0 Ry were fitted up 
to I = 6. The calculations were performed in the local 
density approximation (LDA)iiS and the parametrization 
used was the one of von Barth-Hedin^^. Jomard et a/.— 
showed that generalized-gradient corrections have to be 
included in order to obtain a correct description of the 
stability of the different phases of pure Zr, but as we were 
interested only in the Al-rich part of the phase diagram 
we did not include these gradient corrections. 

Ground state energies at equilibrium Eq, equilibrium 
volumes per atom Vq, and bulk moduli B were obtained 
by fitting the Rose equation of state^i to the energies 
calculated for different volumes around the minimum 

i^(0-i^o(l + ^)exp(-!:^), (1) 

where r is the Wigner-Seitz radius associated to the 
atomic volume V and 5 is related to the bulk modulus B 
through the relation 



For the different compounds, the energies were opti- 
mized with respect to the volume and all other degrees 
of freedom allowed by the symmetry, like the shape of the 
unit cell or some atomic positions. The equilibrium vol- 
umes Vq, bulk moduli B, and formation energies 
relative to the fee phases of both pure Al and Zr are 
presented in table We note that all the formation en- 
ergies are negative, and they thus characterize Al-Zr as 
an ordering system. 

We examined more closely the stability of the phases 
LI2, DO22, and DO23 of AlaZr according to relaxations. 
LI2 being cubic, its energy was optimized with respect 
to atomic volume only, whereas in the tetragonal DO22 
phase optimization was performed additionally with re- 
spect to the c/a ratio and in the tetragonal DO23 phase 
to the c/a ratio and to the atomic displacements Sai and 
Szi {Sa\ and Szi are defined on Fig. 

Our results for AlaZr agree really well with the exper- 
imental ones (Table The equilibrium volumes ob- 
tained for the LI2 and DO23 phases are lower than the 
experimental ones, but this is a known feature of LDA. 
This can be improved by adding gradient corrections: 
Colinet and Pasturel-^ using the generalized gradient ap- 
proximation instead of LDA found a better agreement 
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TABLE I: Equilibrium volume Vo, bulk modulus B, and formation energy E^°^"^ for relaxed Al-Zr compounds calculated with 
LDA. 





Pearson 


Structure 


Vo 


B 






symbol 


type 


(A^/atom) 


(GPa) 


(mRv/atom) 


AI (fee) 


cF4 


Cu 


15.82 


80.78 


0. 


AlaiZr 


cP32 


? 


15.99 


82.56 


-3.04 


AlisZr 


cI32 


? 


16.10 


84.27 


-6.99 


AlsZr 


tI18 


V4Zn5 


16.25 


86.24 


-9.77 


AlrZr (Dl) 


cF32 


CarGe 


16.30 


87.32 


-14.33 


AUZr (Dla) 


tllO 


MoNi4 


16.58 


92.10 


-21.12 


AlaZr (LI2) 


cP4 


CusAu 


16.12 


99.59 


-39.00 


AlaZr (DO22) 


tI8 


AlsTi 


16.60 


99.65 


-39.04 


AkZr fD023) 


tI16 


AlaZr 


16.35 


100.05 


-40.72 


Al2Zr (a) 


hP3 


Cdl2 


18.01 


87.16 


-11.73 


AkZr (IS) 


tI6 


MoSi2 


17.13 


96.40 


-26.19 


AbZr (7) 


0I6 


MoPt2 


17.15 


96.51 


-26.08 


AlZr (Llo) 


tP4 


AuCu 


18.15 


103.33 


-37.07 


AlZr (Lli) 


hR32 


CuPt 


19.04 


93.29 


-16.50 


AlZr (CH40) 


tI8 


NbP 


18.52 


100.48 


-33.56 


AlZr (D4) 


cF32 


7 


18.49 


92.58 


-14.78 


AlZr (Z2) 


tP8 


? 


18.63 


99.70 


-21.03 


Zr2Al (q) 


hP3 


Cdl2 


20.38 


98.10 


-10.72 


Zr2Al (/?) 


tI6 


MoSi2 


19.53 


104.84 


-24.78 


Zr2Al (7) 


0I6 


MoPt2 


19.44 


104.05 


-26.36 


ZrgAl (LI2) 


cP4 


CuaAu 


19.71 


107.67 


-27.11 


ZrgAl (DO22) 


tI8 


AlsTi 


19.88 


105.14 


-23.49 


ZrgAl (DO23) 


tI16 


AlsZr 


19.80 


102.77 


-25.18 


Zr4Al (Dla) 


tllO 


MoNi4 


20.31 


99.85 


-16.30 


ZryAl (Dl) 


cF32 


CarGe 


20.93 


101.66 


-7.96 


Zr (fee) 


cF4 


Cu 


21.70 


98.74 


0. 



"Reference phases are Al(fcc) and Zr(fcc). 



TABLE II: Calculated volumes Vb, c'/a ratios (c' = c/2 for DO22 phase and c' = c/4 for DO23 phase), atomic displacements 
(normalized by c'), and ground state energies relative to the DO23 phase for AlsZr compared to previous calculations and 
experimental data. 





Method 


Vo 

(A^/atom) 


c'/a 


Atomic 
displacements 


AS 
(mRy/atom) 


L12 


Present work 
FP-LMTQ3^ 
VASP2i 
Experiments- 


16.12 

17.4 
17.14 






1.71 
0.64 
2.3 
1.69 


D022 


Present work 
FP-LMTOS 


16.60 


1.141 




1.68 
2.63 




VASB^a 


17.7 


1.141 




1.9 


D023 


Present work 


16.35 


1.087 


(5ai = -0.0021 
<5zr = -0.0273 


0. 




FP-LMTOS 


16.28 


1.09 


(5ai = +0.003 
(5zr = -0.026 


0. 




VASB2^ 


17.5 


1.079 


5ai = +0.0003 
(5zr = -0.0101 


0. 




Experiments— 


17.25 


1.0775 


(5ai = +0.0004 
<5zr = -0.0272 





with experimental data for these equilibrium volumes. 
After relaxing all the degrees of freedom, we see that 
DO23 is the stable phase. As been shown previously by 
Amador et al?'^ using the FP-LMTO technique too, it is 



not enough to consider only the relaxation of the shape 
of the unit cell (c/a ratio) of the phase DO23 to sta- 
bilize it, the atomic displacements 5ai and (5zr allowed 
by the symmetry have to be relaxed too, otherwise LI2 
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still has a lower energy. This was confirmed by Colinet 
and Pastured— with calculations in the pseudopotential 
method, and we observed such a behaviour too. The 
values obtained after relaxation for these displacements 
are close to those measured by neutron diffraction by 
Amador et al?^: the sign of 5a\ is wrong but this rela- 
tive displacement is too close from to be really signif- 
icant. The enthalpy of transformation from the LI2 to 
the DO23 structure was measured by Desch et alm^. The 
experimental value {AH = —1.69 mRy/atom) agrees re- 
ally well with the value obtained from ours calculations 
{AH = —1.72 mRy/atom), which was not the case of 
previous calculations. 

For ZraAl, we found the phase LI2 to have the lowest 
formation energy compared to the two other structures 
we investigated. This is in agreement with the experi- 
mental fact that LI2 is the stable phase of ZrsAl. For 
other compositions, the experimental stable structures 
are not based on an fee underlying lattice, and therefore 
no direct comparison can be made with our calculations. 



III. CLUSTER EXPANSION OF THE 
FORMATION ENERGY 

The FP-LMTO method only allows one to calculate 
the energy of perfectly ordered systems which contain 
a few atoms per unit cell. Disordered or partially or- 
dered systems can be modeled by super-cells, but this 
requires a too large computational time. Moreover, to 
compute the free energy of these systems, one needs to 
calculate the energy of every configuration. This cannot 
be directly done with ab-initio calculations and a cluster 
expansion has to be used. That is why in the following 
we will directly use the FP-LMTO calculations only for 
the perfectly ordered compounds AlsZr in the structures 
LI2 and DO23, and for the solid solution Al-Zr we will 
make a cluster expansion. 



A. Formalism 

Considering a binary crystal of N sites on a rigid lat- 
tice, its configuration can be described through an Ising 
model by the vector <t = {cti, (T2, . . . , cat} where the 
pseudo-spin configuration variable ap equals ±1 if an A 
or B atom occupies the site p. Any structure is then de- 
fined by its density matrix p'', p^{(t) being the probability 
of finding the structure s in the configuration <j. 

To any cluster of n lattice points a — {ii, 12, ■ ■ ■ Tin} is 
associated the following multisite correlation function 



c=Trp^n'^'=^E'^'Hn^- 



(3) 



where the sum has to be performed over the 2^ possible 
configurations of the lattice. 

Clusters related by a translation or a symmetry op- 
eration of the point group of the structure have the 



same correlation functions. Denoting Da the number 
of such equivalent clusters per lattice site, or degeneracy, 
any configurational function can be expanded in the 
fori 



J2b 



dm: 



(4) 



where the sum has to be performed over all non equiv- 
alent clusters and the coefficients fa are independent of 
the structure. 

The usefulness of the expansion Q rests on the fast 
convergence of these coefficients with the size of the clus- 
ters, i.e. with the number of points included in the clus- 
ter as well as the maximal distance between points inside 
the cluster. This allows one to truncate the sum using 
only a finite number of clusters. Knowing the value of 
the function /* for a finite set of structures, the coeffi- 
cients fa can then be obtained using the inverse method 
proposed by Connolly and Williams^S^ i.e. by a matrix 
inversion if the number of structures is the same as the 
number of clusters used in the truncated expansion. Here 
we used more structures than clusters and obtained the 
coefficients by a least square fit. We can thus check the 
convergence of the expansion by its ability to reproduce 

r- 

Rather than doing the fit directly on the configura- 
tional function we did it on the excess function asso- 
ciated which is defined as 



Af = r 



(5) 



where Ci is the point correlation and f^ and f^ are 
the values of the function for a lattice occupied by 
respectively only atoms A {(i = 1) and atoms B (Ci = 
— 1). In the case of the energy, this excess function is 
nothing else but the formation energy. 
Using the expansion 10} we obtained 



Ar 



J2 

(,|a|>2 



o ''1 o 



(6) 

where \a\ denotes the number of points contained in the 
cluster a. 

Applying the Connolly Williams method to the expres- 
sion lO rather than to (0J) allows one to impose easily the 
condition that Af^ should be equal to zero for pure A 
and pure B. We thus obtain the coefficients fa only for 
clusters containing more than one point, the coefficients 
/o and /i of the empty and point clusters being then 
obtained by the relations 



/o 



fA + fB 



2 l±(-l)^Z,„/„, ,7a, 



Q,|a|>2 



h = - E '-^Dafa. (7b) 

Q,|a|>2 
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B. Relaxations 

The volume of a structure, like any other property, 
depends on its configuration. But this volume enters di- 
rectly in the cluster expansion as the coefficients /„ have 
to be calculated for a given volume. As we are gener- 
ally interested in the equilibrium properties, this leads to 
some relaxations. In this study we will consider in two 
different ways these volume relaxations, the globally and 
totally relaxed expansionaSLSS. 

We first can make the cluster expansion explicitly vol- 
ume dependent, writing 



(8) 



where the coefficients fa{V) are obtained by calculating 
the property for different structures at the same vol- 
ume V, the other degrees of freedom (shape of the unit 
cell and positions of atoms) being relaxed, and then by 
using the Connolly- Williams method. Doing such a clus- 
ter expansion for the energy, we can then deduce the equi- 
librium volume associated with a given configuration by 
minimizing with respect to the volume its energy as given 
by expression This is known as the globally relaxed 
scheme and is based on the assumption that the volume 
occupied by every cluster is independent on its configura- 
tion. Such an assumption is questionable in cases where 
there is a significant difference between the atomic vol- 
umes of the constituent elements as in the Al-Zr system. 

Another way to consider relaxations of the volume is to 
calculate the coefficients fa from the equilibrium values 
/*(VJj*) where everything included the volume is allowed 
to relax. The coefficients are then volume independent 
and the values predicted by the expansion are directly the 
ones at equilibrium. Such a treatment is called a totally 
relaxed expansion. This expansion is still rigorous from 
a mathematical point of view since the relaxations are 
themselves functions of the configuration, so the relaxed 
structures will also be. 



C. Results 

1. Total relaxations 

We first made a cluster expansion of the equilibrium 
volume, the bulk modulus, and the formation energy for 
the Al-Zr system on an fee lattice: we are thus doing 
three different totally relaxed expansions. To perform 
the least square fit of the expansion, we used the 26 
structures for which these equilibrium properties were ob- 
tained from our FP-LMTO calculations (Table ID). The 
best agreement between our ab-initio calculations and 
their expansion has been obtained when using 17 clusters: 
the empty cluster {0}, the point cluster {!}, the pairs 
from first to seventh nearest neighbours {2,1}. ..{2,7}, 
seven triangles {3,1}. .. {3,7} presented in figures |21 a- 
g and the tetrahedron of first nearest neighbours {4,1} 



(Fig. Hlh). The values of the coefficients obtained for 
these three totally relaxed expansions are presented in 
table mil and the errors compared to the direct calcula- 
tion for the 25 structures in table lIVl 



(a) {3,1} (b) {3,2} (c) {3,3} (d) {3,4} 



(e) {3,5} 



(f) {3,6} 



(g) {3,7} 



(h) {4,1} 



FIG. 2; Definition of the three and four points clusters on an 
fee lattice used for the expansion. 



TABLE III: Cluster expansion of the equilibrium volume (co- 
efficients Va), bulk modulus (Be), and formation energy (Ja) 
in the total relaxation scheme. 



Cluster 




(A^/atom) 


(GPa) 


(mRy/atom) 


{0} 


1 


18.587 


98.15 


-625.05 


{1} 


1 


-3.230 


-11.12 


419.87 


{2,1} 


6 


0.149 


-1.12 


3.69 


{2,2} 


3 


-0.128 


1.88 


-3.86 


{2,3} 


12 


-0.013 


-0.09 


0.07 


{2,4} 


6 


-0.027 


-0.11 


-0.15 


{2,5} 


12 


-0.037 


0.07 


0.16 


{2,6} 


4 


0.009 


-0.30 


0.93 


{2,7} 


24 


0.014 


-0.17 


0.18 


{3,1} 


8 


0.013 


-0.14 


1.74 


{3,2} 


12 


-0.031 


-0.40 


-0.45 


{3,3} 


24 


-0.001 


-0.14 


-0.55 


{3,4} 


6 


0.023 


0.21 


-0.31 


{3,5} 


24 


0.004 


0.21 


-0.33 


{3,6} 


24 


0.010 


0.14 


0.22 


{3,7} 


24 


0.005 


0.06 


0.54 


{4,1} 


2 


0.030 


-0.55 


0.88 



Looking at the cluster expansion of the formation 
energy, it can be seen that the maximum difference 
between the energy given by the expansion and the 
one directly obtained from the FP-LMTO calculations 
is 4.0 mRy/atom and that the standard deviation is 
1.4 mRy/atom. We did not manage to find a better 
set of clusters producing a smaller error: as we still had 
more structures to fit than clusters, we tried to include 
more clusters like the pair to the eighth nearest neigh- 
bour but this did not improve the difference between our 
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TABLE IV: Deviations for the cluster expansion of the equi- 
librium volume (SVo), the bulk modulus (SB), and the forma- 
tion energy {SE-^°'~"^) in the total relaxations scheme. 

SVo SB SEf°''"^ 
(A^/atom) (GPa) (mRy/atom) 

Ai (fee) a a 0. 

AlaiZr 0.052 0.30 0.42 

AlisZr 0.035 0.83 -1.02 

AlgZr (NbNig) 0.157 -1.05 4.01 

AlrZr (Dl) 0.018 0.31 0.87 

AUZr (Dla) -0.012 -1.62 1.96 

AlsZr (LI2) -0.109 -0.02 -0.21 

AlsZr (DO22) -0.029 0.63 -0.58 

AlsZr (DO23) -0.079 0.73 -2.10 

AhZi (a) -0.017 0.10 -0.36 

AbZr (/?) -0.028 0.52 -1.11 

Al2Zr (7) -0.041 0.60 -1.38 

AlZr (Llo) 0.219 -1.18 1.96 

AlZr (Lli) 0.323 -0.86 0.60 

AlZr (CH40) 0.077 -0.11 0.56 

AlZr (D4) -0.348 0.64 -1.20 

AlZr (Z2) 0. 0.41 0.49 

Zr2Al (a) 0.011 0.03 0.16 

Zr2Al 1(3) -0.007 0.46 -0.71 

ZraAl (7) 0.002 0.49 -0.59 

ZraAl (LI2) -0.095 1.21 -0.32 

ZraAl (DO22) -0.038 2.01 -1.29 

ZraAl (DO23) -0.061 -2.03 -0.69 

Zr4Al (Dla) -0.012 -1.62 1.96 

ZryAl (Dl) 0.085 0.32 1.78 

Zr (fee) 0. 0. 0. 

Standard deviation 0.116 0.91 1.33 



FP-LMTO calculations and their expansion. Such an er- 
ror does not allow one to reproduce the relative stability 
between different ordered compounds at a same concen- 
tration, for instance between the phases LI2, DO22, and 
DO23 of AlaZr. But as we are interested in using the 
cluster expansion only for the solid solution Al-Zr, this is 
not a problem: for perfectly ordered compounds we can 
directly use the results of our ab-initio calculations. 

For the equilibrium volume, we can compare the accu- 
racy of the cluster expansion with the one of the Vegard's 
law which assumes a linear relation between the atomic 
volume and the concentration. The standard deviation of 
the Vegard's law is 0.427 A/atom. For none of the con- 
sidered structures such an important error occurs, and 
we have thus obtained an important improvement by not 
considering only the empty and point cluster as one does 
in the Vegard's law. 

For the bulk modulus, the accuracy of our FP-LMTO 
being of the order of 1 GPa, here too we can consider the 
convergence of the cluster expansion to be good. 



2. Global relaxations 

Using the same sets of clusters and structures, we ex- 
panded the ground state energy in 21 different volumes 
between 14 and 24 A'^/atom. For each structure, these 
21 expansions gave the ground state energy of the re- 
laxed structures at the corresponding fixed volume, and 
we then used these results to obtain the volume, the bulk 
modulus, and the ground state energy at equilibrium by 
fitting the Rose equation of state2i. 



TABLE V: Deviations for the cluster expansion of the equilib- 
rium volume (SVo), the bulk modulus (SB), and the formation 
energy [5E^°^"^) obtained in the global relaxations scheme. 





SVo 


SB 






(A^/atom) 


(GPa) 


(mRy/atom) 


Maximal deviation 


-0.351 


-2.26 


5.15 


Standard deviation 


0.120 


0.94 


1.44 


The maximal and 


standard 


deviations 


between the 



properties deduced from the expansion and the ones 
directly obtained from the FP-LMTO calculations are 
shown in tabled They are close to what we previously 
obtained in the total relaxation scheme. Actually, we 
did not get any important difference between the results 
obtained according to the way volume relaxations are 
considered. For each structure the deviation is quite the 
same in both cases, this being true for the formation en- 
ergy as well as for the equilibrium volume and the bulk 
modulus. Such a result could not have been easily pre- 
dicted as the size difference between Al and Zr is quite 
important: the atomic volumes given by our calculations 
for the fee structures of Al and Zr are respectively 15.82 
and 21.70 A^ (Table H)). This proves that the globally 
and locally relaxed expansions are equivalent. 

As the totally relaxed expansion only gives us one set 
of coefficients for the whole range of volumes, it is more 
convenient and we will use this expansion in the follow- 
ing. 

IV. FINITE TEMPERATURE PROPERTIES 

At finite temperature, the vibrational and electronic 
contributions as well as the configurational entropy have 
to be included in the description of the system. Consid- 
ering two different time scales, a slow one for the con- 
figurational effects and a much faster for vibrations and 
electronic excitations^, we define a vibrational and elec- 
tronic free energy, F"*''(cr) and F'^'(<t), both depending 
on the configuration. Using the variational principle, the 
free energy is obtained by minimizing the functional 

F[p] = {Eo) + {F^^') + {F'^') + ksTilnp), (9) 

where ks is the Boltzmann constant and p the density 
matrix. 



Al (fee) 


0. 


0. 


AlaiZr 


0.052 


0.30 


AlisZr 


0.035 


0.83 


AlgZr (NbJNig) 


0.157 


— 1.05 


A 1 T~I / 1 \ 

AlrZr (Dl) 


0.018 


0.31 


AUZr (Dla) 


—0.012 


— 1.62 


AlsZr (LI2) 


—0.109 


-0.02 


AlsZr (DO22) 


—0.029 


0.63 


AlsZr (DO23) 


—0.079 


0.73 


AbZr (a) 


—0.017 


0.10 


AbZr (/?) 


—0.028 


0.52 


AbZr (7) 


—0.041 


0.60 


AlZr (Llo) 


0.219 


— 1.18 


AlZr (Lli) 


0.323 


—0.86 


AlZr (CH40) 


0.077 


-0.11 


AlZr (D4) 


-0.348 


0.64 


AlZr (Z2) 


0. 


0.41 


Zr2Al (a) 


0.011 


0.03 


Zr2Al (/3) 


-0.007 


0.46 


ZraAl (7) 


0.002 


0.49 


ZraAl (LI2) 


-0.095 


1.21 


ZraAl (DO22) 


-0.038 


2.01 


ZraAl (DO23) 


-0.061 


-2.03 


Zr4Al (Dla) 


-0.012 


-1.62 


ZryAl (Dl) 


0.085 


0.32 


Zr (fee) 


0. 


0. 


Standard deviation 


0.116 


0.91 
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The cluster expansion of the formation energy at T = 
K gives us an expression for the cohesive part of the 
functional of Eq. 0. We do not have to take into ac- 
count any variation of the lattice parameter with temper- 
ature as we chose to work in the harmonic approximation: 
Ozolins and Astci showed on the solubility limit of Sc in 
Al that there was only a small improvement when going 
from the harmonic to the quasiharmonic approximation. 
Similar expressions have to be found for the electronic 
and vibrational parts of the expression ©. The mini- 
mization of F[p\ with respect to p will then be done in 
the Bragg- Williams approximation. 

A. Electronic free energy 



Zr (fee), AlgZr (LI2), and AlgZr (DO23) (Fig. EJ. In the 
range of temperature of interest, i.e. below 1000 K, this 
electronic contribution is smaller than 1 mRy/atom, and 
so is the excess free energy associated. This is the same 
range of order as the accuracy of the cluster expansion 
of the formation energy. We thus chose to neglect this 
contribution to the free energy. 

B. Vibrational free energy 

We studied the vibrational effects in the harmonic ap- 
proximation, comparing the ability of the Debye model 
with a phonon calculation to predict the thermodynamic 
properties. 



At a temperature of K, all electronic states of en- 
ergy below the Fermi level ef are occupied, whereas the 
ones above are empty. At finite temperature, the elec- 
trons close to the Fermi levels can be promoted to states 
of higher energies according to the Fermi-Dirac distribu- 
tion f(e,T). The electronic excitations induce a change 
of the charge density and thus of the effective potential of 
the one electron Hamiltonian. This leads the electronic 
density of states (DOS) n(e) to be temperature depen- 
dent. But the changes induced on the total energy and 
on the entropy by this temperature dependence of the 
electronic DOS are small^. We thus assumed the elec- 
tronic DOS to be temperature independent and equal to 
the one obtained at T = 0. The energy change l^E"^ (T) 
and the entropy S'^^{T) due to electronic excitations are 
then 

/oo 
e7i(e)[/(e,T)--/(e,0)]de, (10a) 
-00 

/OO 
n(e){/(e,r)ln[/(e,T)] 
-00 

+ [l-/(e,r)]ln[l-/(e,T)]}de. (10b) 



(d) 



Ja) 



-lb) 



Al(fcc) (a) 

Zrffcc) b) 

Al3Zr(Ll2) (cJ 

Al3/r(D023) (d) 



300 400 500 600 700 BOO 900 1000 
T(K) 



FIG. 3: Electronic free energy, F-^' = AE"' - TS'''. 



We calculated the electronic contribution to the free 
energy, F"^ = A£"=' - TS''\ for the structures Al (fee), 



1. Phonon calculation 

A calculation of the phonon DOS n{uj) allows one to 
compute the vibrational free energy. For temperatures 
higher than 300 K, it is enough to consider only its high 
temperature expression 



F 



'vih 



+0 



/>oo 

— 31n(A;Br)+ / ln(fiu))n(a;)da 

^0 



T 



(11) 



Phonon DOS were calculated for Al (fee), Zr (fee), 
AlaZr (LI2), and AlaZr (DO23) in the linear response 
theory frameworkSS. We used energy independent MTO 
as a basis for representing the first order correction to 
the one electron wave functions in the implementation 
developed by Savrasov'^"''^^. These calculations were per- 
formed in the LDA using the parametrization of Moruzzi- 
Janak- Williams^. The radii of the muffin-tin spheres 
were taken equal to the ones of the band structure cal- 
culation. For valence states, the basis used was the same 
whereas the 4s and 4p states of Zr were treated in two 
different panels with the respective kinetic energies 
2.7 and 1.1 mRy. For fee structures, phonon frequencies 
were calculated on a grid of 8 x 8 x 8 wave vectors q lead- 
ing to 29 points in the irreducible Brillouin zone (IBZ), 
for LI2 a grid of 5 x 5 x 5 leading to 10 points in the IBZ 
was used, and for DO23 a gi'id of 4 x 4 x 4 wave vectors 
leading to 13 points. 

The calculated phonons dispersion for Al fee is com- 
pared in figure ^ to the measurements of Stedmam et 
^^j33j34 -fpj. ^]-^j.gg different high-symmetry directions. We 
see that our calculation overestimates phonon frequency. 
Other phonon calculations for Al fcoSi^Si2i using the lin- 
ear response theory too obtained a better agreement with 
experimental data. They all used a plane waves basis in 
the pseudopotential framework, but the use of an energy 
independent MTO as a basis does not seem to be the 
reason of the discrepancy with experimental data in our 
case, as Savrasov showed for Nb^Sas well as for NbC and 
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FIG. 4: Calculated phonons dispersion for Al fee in solid line 
compared to experimental data (Refi— i— ) denoted by crosses. 



Si^ that this basis was well-suited to obtain phonon dis- 
persion. 

The phonon DOS obtained from these calculations for 
Al (fee), Zr (fee), AlgZr (LI2), and AlgZr (DO23) are 
presented in Fig. \S\ For Al (fee), we compared our cal- 
culated phonon DOS with experimental ones. Experi- 
mental DOS^^i22i^ were obtained by means of a Born- 
von Karman model. Force constants were fitted up to 
the eighth nearest neighbours in order to reproduce the 
phonon measurements in high-symmetry directions of 
Stedmam et alm^^, the Born-von Karman model be- 
ing used then to compute the frequency distribution. We 
can see on the phonon DOS too that our calculated fre- 
quencies are slightly too high. Nevertheless, the shape of 
the frequency distribution is correct. 



2. Debye model 

The Debye model assumes a linear dispersion between 
the phonon frequency and its wave vector. This leads to 
the following high temperature expression of the vibra- 
tional free energy 



= knT 



-1 



(12) 



where the Debye temperature 0£) is obtained from the 
elastic constants of the structuresiS-. 

The elastic constants were obtained by means of FP- 
LMTO calculations using the same set of parameters as 
for the formation energy calculations. The unit cell of 
the crystal was deformed around its equilibrium posi- 
tion in order to obtain the second derivative of the en- 
ergy at its minimum which can be then related to the 
elastic tensori^ii^. During this deformation, no relax- 
ation was allowed. For the DO23 structure, the c/a ratio 
and the position 5ai and dzr of the atoms were frozen at 
their equilibrium value. For some of the deformations, we 
checked that these relaxations did not change much the 



values of elastic constants. Moreover, as we are lowering 
the symmetry of the structure by deforming it, some new 
degrees of freedom can appear, but we did not consider 
either these ones. 

The elastic constants calculated with the FP-LMTO 
are compared to the experimental ones in table rVTl The 
discrepancy between the calculated and experimental 
constants is in the order of 10%. This leads to some dif- 
ferences between the Debye temperatures obtained from 
these calculated constants and the ones obtained from 
the experimental constants, but the relative positions of 
these temperatures are correctly predicted. 

In table IVIII we show Debye temperatures obtained 
from a calculation of the elastic tensor for cubic struc- 
tures fee, Dl, and LI2 of the Al-Zr system. The structure 
D4 of AlZr is cubic too, but this phase was found to be 
mechanically unstable through a Bain deformation path 
and cannot be used to calculate a Debye temperature. 



3. Comparison for ordered compounds 

As we calculated the phonon spectrum for AlsZr for 
the stable structure DO23 and the metastable one LI2, we 
were able to compare the excess vibrational free energy 
^pvib obtained from the phonon DOS and the Debye 
model, the reference phases being Al (fee) and Zr (fee) 
(high temperature expressions are given in table IVIII|I . 
We thus see that the Debye model makes an important 
error in predicting this excess free energy as it overes- 
timates it by a factor ^ 2. This error comes from the 
inability of the Debye model to reproduce the phonon 
DOS as shown by figure |31 Moreover the phonon calcu- 
lation shows that the two considered structures of AlaZr 
should have the same vibrational free energy which is 
not correctly predicted by the Debye model. This error 
of the Debye model would lead to a stabilization of the 
phase LI2 at high temperatures (T >905 K) which is not 
true experimentally. In order to correctly describe the 
relative stability of these two ordered phases of AlaZr we 
cannot use the Debye model, but we need the phonon 
calculation. 



4- Cluster expansion for the disordered phase 

For the vibrational free energy of the disordered phase, 
we made a cluster expansion of the vibrational free en- 
ergies of several ordered structures. As the Debye model 
only requires the calculation of the elastic tensor, which is 
much more faster than a calculation of the whole phonon 
spectrum, we used it to calculate the vibrational free en- 
ergy of these ordered compounds (Debye temperatures 
in Table IVIII) . By doing so we saw previously that we 
overestimates AF'"'^^, but a calculation of the phonon 
spectrum is not conceivable for a number of structures 
large enough to fit the cluster expansion. We have then 
to accept such an error. 
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(e) AlgZr (LI2) 



(d) AlaZr (DO23) 



FIG. 5: Phonon density of state. 



TABLE VI; Elastic constant dj (in GPa) calculated with the FP-LMTO compared to experimental values for Al (fee), AlsZr 
(D023), and Zr (hep), and Debye temperature 9d associated. Debye temperatures obtained by calorimetric measurements of 
the specific heat, when available, are given in brackets. 







Cii 


C33 


C12 


Cl3 


C44 






Bd (K) 


Al (fee) 


FP-LMTO 


10L5 




70.4 




31.7 




385 






exp4S^ 


114.3 




61.9 




31.6 




431 


(428)^ 


AlgZr (D023) 


FP-LMTO 


215.3 


228.2 


54.1 


33.3 


103.2 


123.5 


616 






exp^i'' 


208.8 


208.3 


70.5 


49.1 


87.2 


102.2 


575 




Zr (hep) 


FP-LMTO 


153.1 


171.2 


63.4 


76.5 


22.4 


44.9 


262 






exp.«'*« 


155.4 


172.5 


67.2 


64.6 


36.3 


44.1 


299 


(310)^ 



"measured at room temperature 



Looking at the high temperature expression of the vi- 
brational free energy given by the Debye model (Eq. I12|l . 
we can make the following cluster expansion 



which allows us to write the vibrational free energy as 



F 



•vih 



^DaJaCa - 3 In T 



(14) 



31n6l 



D 



(13) 



By doing so, the temperature dependence of the free en- 
ergy is really simple and we do not have to make a cluster 
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TABLE VII; Elastic constants dj (in GPa) for Al-Zr com- 
pounds of cubic symmetry and Debye temperature 6d asso- 
ciated. 





Cii 


C12 


C44 


en (K) 


Al (fee) 


101.5 


70.4 


31.7 


385 


AlyZr (Dl) 


136.5 


62.7 


45.8 


449 


AlgZr (LI2) 


187.3 


55.7 


95.1 


557 


ZrgAl (LI2) 


163.8 


79.3 


86.5 


388 


ZryAl (Dl) 


136.3 


84.4 


56.6 


300 


Zr (fee) 


121.4 


87.1 


45.7 


249 



TABLE VIII: Comparison of the high temperature expres- 
sions of the vibrational free energy obtained with the phonon 
calculation and the Debye model. 



AlgZr (LI2) 


Phonons 




= 0.85A;Br^ 


hO(^) 




Debye 




= 1.44fcflr- 


hO(^) 


AlgZr (DO23) 


Phonons 




= 0.85fcsr- 


hO(^) 




Debye 




= 1.74fcsr^ 


hO(i) 



expansion of the free energy at every temperature. 

We only used four clusters in the truncated expansion: 
the empty cluster {0}, the point cluster {1}, the pair 
{2,1} of first nearest neighbours, and the triangle {3,1} 
of first nearest neighbours. The eight structures of table 
IVIII were used to fit the coefficients of the expansion. The 
results of this expansion are presented in table llXt a) and 
the deviations in table llXt bl. Although only few clusters 
were used in this expansion, the convergence is really 
good. 



TABLE IX: Cluster expansion of the function f = 3 In So - 1 
for the vibrational free energy. 











r 


5r 


Cluster 




J a 


AI (fee) 


16.86 


0. 


{0} 


1 


17.385 


AlrZr (Dl) 


17.32 


-0.08 


{1} 


1 


0.874 


AlsZr (LI2) 


17.97 


0.04 


{2,1} 


6 


-0.197 


ZrsAl (LI2) 


16.88 


0.04 


{3,1} 


8 


-0.027 


ZryAl (Dl) 


16.11 


-0.08 








Zr (fee) 


15.55 


0. 



(a) Coefficients of the 

expansion. (b) Deviation of the 

expansion. 



C. Bragg- Williams approximation 

We thus obtained an expression for the different parts 
of the free energy functional F[p\ of expression ©: the 
cohesive part is given by the cluster expansion of the 
FP-LMTO calculations (coefhcients in table IIII|) . the vi- 
brational energy by the expression H14|) with the coeffi- 
cients of table HxT a). and the electronic contribution can 



be neglected. The functional F[p\ is minimized in the 
Bragg- Williams approximation. This assumes that there 
is no short range order and that the correlation functions 
can be factorized over the mean values of the pseudo spin 
variable (ct„) for the lattice sites contained in the cluster, 

Cc. = ([{a,) = X{{o^). (15) 

The Bragg- Williams approximation thus assumes that 
the lattice sites interact only through their mean oc- 
cupancy and neglects all correlations between different 
sites. This can be improved by using the Cluster Varia- 
tion Method (CVM)^, but in the case of a low solubility, 
no really important improvement is expected when going 
from the Bragg- Williams approximation to the CVM. 
Moreover, the computational time necessary to obtain 
the free energy by means of the CVM increases a lot 
with the size of the maximal cluster. As Zr has a really 
low solubility in Al (fee) and as the long range interac- 
tions of the cluster expansion of the formation energy 
requires a too large cluster, we chose to work with the 
Bragg- Williams approximation. 

Within the Bragg- Williams approximation, the con- 
figurational entropy has the following expression for a 
binary compound 

s[p\ = -fcs5I(i + K))in(i + K)) 

n 

+ (l-(a„))ln(l-(a„)). (16) 

1. Disordered phase 

For a disordered state, all lattice sites are equivalent 
by symmetry. They have thus the same point correlation 
Ci — 2x — 1, where x is the Zr atomic concentration. 
Consequently any correlation function can be written in 
terms of the point correlation: 

Ca=Cl"'. (17) 

The cluster expansion of the function using the ex- 
pression © of the excess function A/", can then be ex- 
pressed as a function of the point correlation, or equally 
as a function of the concentration. This leads to an ex- 
pression similar to the way the internal energy of a solid 
solution is written in a Redlich-Kister model which is of 
common use in the Calphad method^ 

r = xf^ + (1 - x)/^ + x{l -x)Y, Ln{2x - 1)", (18) 

n>0 

where the coefficients L„ are obtained from the coeffi- 
cients Ja by the relations 

i>l , a 

- \a\=n+2i 
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Using the expression H16|l for the entropy, we obtain 
for the free energy of the disordered fee sohd solution 
Al(i_x)Zra; 

F{x) = {1 - xp^^^^"" + xU^'''^'"' 

+kBT[xlnx + {1 -x) ln(l ~ x)] 

+a;(l-a:)^L„(2a;-l)". (20) 

n>0 

The Redhch-Kister coefficients are obtained from the 
cluster expansion of the formation energy and the cluster 
expansion of the vibrational free energy calculated in the 
Debye approximation. 

Lo = -89.09 + 29.9 x 10~^T mRy/atom 

Li = -14.30 + 5.47 X lO^^T mRy/atom (21) 

L2 = —7.03 mRy/atom 

For a dilute solution {x ^ 1), the expression f^^i is 
equivalent to the free energy of a regular solution, the 
excess free energy being then x{l — x)n = x{l — x){Lo + 
Li +L2)- We compare in table|X]the value of obtained 
from our calculations to the ones obtained by a fit of the 
phase diagram through a Calphad approach^i*^. 

TABLE X: Parameter Q, (in mRy/atom) of the excess free 
energy for an fee regular solid solution Al-Zr deduced from 
ab-initio calculations and compared to values obtained by a 
fit of the experimental phase diagram. 



Present work 


n 




-110.42 + 35.37 X 10 


— Zrj-i 


Saundera^ 






-87.60 + 22.85 x 10" 




Murray et al}- 






-85.08 + 31.01 X 10" 





2. Line compounds 

AlsZr in the DO23 or LI2 structures can be considered 
as a line compound, i.e. perfectly ordered: both struc- 
tures are composed of interpenetrating sublattices of pure 
Al and pure Zr. The configurational entropy of such line 
compounds can be neglected and these structures only 
exist for a concentration x = 1/4. We checked with a 
calculation using the previous cluster expansions of the 
formation and vibrational energies that this assumption 
was correct in the range of temperature of interest when 
looking at equilibrium with the solid solution. The free 
energy of these compounds is then simply given by 
Q 1 

pAhZr _ _jjAlJcc _|_ _jjZr,fcc ^ ^-yAl-iZr (22) 

4 4 ' V ; 

where Af/^^'^Zr ( 

in mRy/atom) is obtained from our pre- 
vious calculations of the formation energy (Tabled and 
of the excess vibrational free energy calculated from the 
phonon DOS (Table IVnH) . 

^jjAi:,Zr,Li2 ^ -39.00 + 5.38 X lO^^T, (23a) 

^^Ai,Zr.D023 = -40.72 + 5.38 x lO^^T. (23b) 



D. Solubility limit of Zr in Al (fee) 

Using the previous expressions for the free energies of 
the disordered phase and the line compounds AlaZr, we 
obtained the solubility limit of Zr in Al (fee), both in 
the stable phase diagram when considering the structure 
DO23 foi' AlsZr and in the metastable one when consid- 
ering the structure LI2. As we are in the case of a dilute 
solid solution, the solubility limit of Zr in Al (fee) is an 
analytic function of the temperature^, 

. = exp(^ k^J^)- (24) 

The solubility we obtained is too low: at the 
melting temperature of the perictectic it is equal to 
0.0016 at.% Zr, whereas the one deduced from exper- 
imental data is 0.08 at.% Zr-S^, When comparing the 
variation with respect to l/T of Inx with experimental 
measurements, we obtain a straight line having the same 
slope as Fink data^ (c/. Fig. 01. This shows that our 
calculations provide an approximation of the enthalpy 
difference between the solid solution and the DO23 or- 
dered compound which is consistent with Finks data, and 
that the discrepancy on the solubility limit only arises 
from an error on the estimation of the entropy differ- 
ence. Computing the solubility limit of Sc in Al, Asta et 
al.^ reached the same conclusion that ab-initio calcula- 
tions correctly predicted the enthalpy difference between 
the ordered compound and the solid solution when com- 
pared to experimental data. In our case, the error on 
the entropic terms may come from an overestimation of 
the vibrational free energy of the disordered phase due 
to the use of the Debye model for this phase. As for the 
structures DO23 and LI2 of Al3Zr, the Debye model over- 
estimates the excess vibrational free energy by a factor 
~ 2 (Table IVIII|I , we can think that we get an error of 
the same range for the solid solution. 
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FIG. 6: Fit of the entropy of the solid solution so as to repro- 
duce Fink experimental data. 

We correct the entropic part, leaving unchanged the 
enthalpic part, of the parameter f2 defining the excess 
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free energy of the solid solution so as to obtain a perfect 
agreement with Fink data (c/. Fig. 1^, and we obtain 

n = -110.42 + 10.07 X lO^^T mRy/atom. (25) 

We thus get a stable solubility limit that is consistent 
with Fink measurements, and we are now able to predict 
the metastable limit using the expression (|25|l to evaluate 
the excess free energy of the solid solution. As the struc- 
tures DO23 and LI2 of AlsZr have the same vibrational 
free energy, the difference of solubility limit is only due 
to the difference of ground state energies of these two 
phases. At the melting temperature of the peritectic, 
we obtain a maximal metastable solubility limit equal to 
0.275 at.% Zr. 
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FIG. 7: Calculated stable and metastable solubility limits of 
Zr in Al compared to experimental data&2*— *— 

This is to compare to the results obtained by a Cal- 
phad method. Murray et alm^ modelled the metastable 
phase of AlaZr as a line compound. They assumed 
that only an enthalpy term, and no entropy term, con- 
tributes to the free energy difference between the stable 
and metastable phases. This was done to guarantee that 
LI2 does not become stable at high temperature. More- 
over there is no experimental data that allows to esti- 
mate the entropy of the LI2 phase. Our calculation of 
the vibrational free energy shows that such an approx- 
imation was correct. The enthalpy difference between 
the two compounds was assumed to arise from the co- 
herency of the phase LI2 with the matrix. They calcu- 
lated from the elastic properties of Al and an estimate 
of the composition dependence of the lattice parameter 
an elastic energy of 1.52 mRy/atom. This estimation is 
quite close to our calculation {AH — 1.72 mRy/atom) 
as well as to the experimental measurement of Desch et 
alm^ {AH = 1.69 mRy/atom). They thus obtained a 
solubility limit that is higher in the metastable phase di- 
agram than in the stable one and their prevision is really 
close to our result: they predicted a maximal metastable 
solubility limit equal to 0.21 at.%. 

In another Calphad study, Saunders-i used the Gibbs 
energy for the disordered Al (fee) solution, as derived 



from the stable equilibrium diagram, to construct the 
Gibbs energy of the ordered LI2 phase in the Bragg- 
Williams approximation. He found a higher solubility 
limit for Zr in the metastable phase diagram than Murray 
et al. as he predicted a metastable solvus composition of 
0.3 at.% Zr at the melting temperature of the peritectic. 

Our study thus allows one to estimate the free energy 
difference between the stable and metastable phases of 
AlaZr, quantity which is not available experimentally and 
has to be guessed in these Calphad studies. One thus 
sees how it is possible to improve the thermodynamic 
database available to Calphad methods. 



V. CONCLUSION 

The equation of state for several compounds in the 
Al-Zr system has been computed using the full potential 
linear muffin tin orbital method (FP-LMTO). These ab- 
initio calculations correctly predict the stability of the 
phase DO23 for AlaZr if we consider the cell internal re- 
laxations. 

We made a cluster expansion of the results of ab-initio 
calculations to predict the formation energy of any com- 
pound in the Al-Zr system based on an fee underlying lat- 
tice. We showed that despite the size difference between 
Al and Zr a totally or globally relaxed expansion for the 
volume leads to the same result: there is no difference if 
we use the cluster expansion to predict the formation en- 
ergy at the equilibrium volume of each structure or at a 
fixed volume, the energy being then minimized according 
to the volume. 

For finite temperature calculations, we showed that 
the electronic excitations can be neglected. The vibra- 
tional energy was studied in the harmonic model, using 
different levels of approximation: the Debye model was 
compared to results obtained from a calculation of the 
phonon spectrum for AlaZr in the structures DO23 and 
LI2, and it was found that the use of the Debye model 
leads to an overestimation of the vibrational free energy. 
So we preferred to use the results from the phonon spec- 
trum to calculate the vibrational free energy of ordered 
compounds. For the disordered phase, we chose to make 
a cluster expansion of the vibrational free energy. It was 
only possible with the Debye model as this requires less 
computational time. 

We were able to calculate the solubility limit of Zr in 
Al (fee) in the the Bragg- Williams approximation. The 
solubility limit obtained is too low compared to experi- 
mental data. We showed that this discrepancy is due to 
an error in the estimation of the entropy in our thermo- 
dynamic model. This may be due to an overestimation 
of the vibrational free energy of the disordered phase 
because of the use of the Debye model for this phase. 
Correcting the vibrational entropy of the solid solution 
so as to fit the experimental measurements of Fink, we 
were able to predict the metastable solubility limit which 
lies between the estimation of Murray and the one of 
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Saunders, both obtained by a Calphad method. We thus 
showed how first principles calculations can lead to an 
estimation of the phase diagram. This approach has the 
advantage of not requiring any experimental input, and 
consequently this is not a problem to predict stability of 
metastable phases. 
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